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A METHOD FOR COMPUTING VERTICAL-PLANE COVERAGE 
DIAGRAMS FOR FREQUENCY AGILE PULSE RADAR SYSTEMS 


INTRODUCTION 


In a comprehensive report on maximum radar range computations, Blake ri] 
provides all the necessary tools for obtaining coverage plots (i.e., range- 
height-angle charts) of monochromatic pulse radar systems. This paper provides 
a method for extending the range computations applicable to the monochromatic 
case to pulse radar systems featuring frequency agility. The number of dif- 
ferent frequencies employed, as well as the number of consecutive pulses to 
be transmitted at the same frequency, are both arbitrary. 


MATHEMATICAL ANALYSIS 


The detection scheme assumed in the analysis which follows utilizes a 
square-law detector followed by a linear integrator, as shown in Fig. 1. The 
symbols Dj, Xo, and Y, denote the SNR (i.e., signal-to-noise power ratio) at 
the designated points of the detector-integrator complex. 


/ Based on approximate closed-form expressions derived by Barton f2] and 
Fe later simplified by Cann Esi; one may relate the input and output SNRs using 
the equation 


(1) 
Dot Zed 


The integrated output SNR, Yo may now be written in the form 


N N 2D... 
° 
Yy= © X= 5 fe (2) 
; J=2 J=2 By + 23 
le y where N denotes the number of integrated pulses and D.. represents the SNR at 
| x the detector input due to the jt pulse illuminating a target located at R,..; 


it is assumed that Ryox represents the maximum detection range associated with 
the total energy delivered on the target by the N transmitted pulses. 


Since our objective is to relate the integrated SNR (i.e., Y,) toa 
maximum detection range, Rmax, uSing the classical radar equation, an average 
SNR is defined at the detector output in the form 


Note: Manuscript submitted March 4, 1977. 
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(3) 


1 j=l Dyj+2.3 


The SNR given by Eq. (3) may be thought of as the SNR necessary to achieve a 
specified probability of detection, Py, and false alarm, Pr,, for a specified 
number of pulses and a postulated type of target; normally, it is the corres- 


a 


ponding SNR, D,, at the input of the square~law detector rather than .. that 
one associates with Py and Pgg.- Given, however, D,, the corresponding value 
of X, is obtained using Eq. (1). For example, given N, Py, and Pras and a 

target described by Swerling's Case 2, one may utilize published results {1] 


to determine Bo and subsequently X,. 


Once x, is obtained, the problem becomes one of relating i. to a corres- 
ponding maximum detection range, Ryax, in a way which takes advantage of the 
computational aids developed by Blake [1,4!. This is accomplished by writing 
the radar equation in the forn, 


2 4 2 
PA oF, P,Ao = B 


ar 28 J 4 
max anh; Do Pnoise Rnax 


(4) 


2 . 
vs Phoise’® 


where Dj; represents the SNR at the detector input due to the jth pulse return 
from a target positioned at Rmax; Dy is the (averaged) minimum SNR based on N 
pulses; F, denotes the pattern-propagation factor; and the remaining symbols 
in Eq. (4) represent well-known radar parameters [5]. The bracketed portion 
of Eq. (4) represents the free-space range, S50," corresponding to N pulses, 
but as if they all had been transmitted at frequency f; (i.e., wavelength Aj) 
and each resulted in a SNR, D,, at the detector input. In view of the above 
interpretation, one may rewrite Eq. (4) in the forn,- 


Eg. \* 
=f xJ 
Do 3 = D, (= ) (5) 


where Ryj = Rits 


Since 


jts ' (6) 


= 3 
Anh, Do Phoise 


varying only the frequency results in a relationship of the form 


dy f * 
a oT - aT at 
Rjts = Riss ae Bits \, (7) 
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which provides the free-space range at frequency f; in terms of the free-space 
range at another frequency, fj, while all other parameters in Eq. (6) remain 


a Se ay 


fixed. 
In some instances the maximum free-space range, R. » based on the 
p number of pulses at each frequency, is known. In this case, there is a Boj 
j given by 
P Avon ,* 
D = ° (8) 
MD ah? P R ‘ 
j ‘noise “j max 


The free-space range at frequency f; is then given by 


ae 
4 _ 4 oj 
Ries a a max ( Bo ) (9) 


Substituting next Eqs. (3) and (5) into Eq. (1), that is, into expression 


2D 
a ° 
X, = (10) 
5 1 2.3/5. 
4 6, 
one obtains 
i (= i 
: N R 
: 1 i" = r ~ (11) 
i= max\4 
1+@& j=l 1 oe 
Rj 
where %= (2,3/D,). 
ig PLOTTING TECHNIQUES 
r Solution of Eq. (11) for R,,, provides the information necessary to 
~ produce the so-called radar vertical-plane coverage diagrams. 


Blake [4] has developed a computer program for presenting the radar 
interference lobing phenomena on range-height-angle plots. This program was 
named LOBEPLOT and was written in Fortran for use on a CDC-3800 computer. 


Relatively simple modifications to LOBEPLOT have been made which have 
resulted in the solution of Eq. (11) and enabled the computer (Calcomp) plotting 
of the lobing charts for frequency agile radars. The original modifications 
were made, and the program was debugged on a CDC=3800 computer at NRL. More 
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recently, the program has been changed to allow running on the new Texas 
Instruments Advanced Scientific Computer (ASC) at NRL. This new version has 
been named LOBMUF for “lobes, multi-frequency." The complete program is 
listed in Appendix A; the input data card formats are given in Appendix B. 


In LOBMUF, Pg, Prag, the total number of pulses integrated, the Swerling 
fluctuation case, and the number of radars or radar frequencies are specified 
on input data cards. In addition, for each radar or radar frequency, the 
free-space range, beamwidth, number of pulses at that frequency, sidelcbe 
level, and a parameter FREF are required. FREF is zero if the free-space 
range was calculated at each frequency as indicated in connection with Eqs. 
(8) and (9); otherwise, FREF is the midband or average frequency and the free- 
space range corresponding to this reference frequency. Other inputs to the program 
are identical to those for LOBEPLOT. They concern the dimensions of the range- 
height chart, polarization radiated, antenna height and tilt, etc. 


Using subroutines from another of Blake's computer programs [5], RGCALC, 
the signal-to-noise ratio, Bo, and the parameter, Q, are computed, If FREF = 0, 
Boj for each frequency is calculated. Next, the pattern propagation factor, 
Fy, is calculated exactly as in Blake's LOBEPLOT. Finally, a simple search 
routine is used to converge on the value of Rmax in Eq. (11). The search is 
terminated when the two sides of Eq. (11) differ by less than 0.01. 


Figure 2 gives a typical coverage diagram generated by the computer program 
LOBMUF. The inputs involved a hypothetical radar radiating two pulses at each 
of 21 frequencies from 1350 MHz toc 1650 MHz in 15 MHz steps. 1500 MHz was 
chosen as the reference frequency, and the free-space range was chosen to be 
100 n. mi. 


Figure 3 shows the lobing plot for a single frequency radar at 1500 MHz. 
A comparison of Figs. 2 and 3 shows that the incoherent integration in a 
frequency agile radar does much to fill in the nulls of the lobing pattern and 
give a solid elevation coverage. 


SUMMARY 


A computer program has been developed to plot range-height-angle lobing 
charts for frequency agile radars. This program is quite flexible. It may 
be used to plot the lobing charts for several radars or a single radar including 
several frequencies. This program, LOBMUF, builds on Blake's range-height 
charts and lobing plots. The program assumes incoherent integration and a 
closed form expression for a square-law detector. 
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Fig. 3 — Lobing chart for single frequency radar at 1500 MHz with free-space range of 
100 n. mi. A wave height of 4 feet and vertical beamwidth of 4.0° were assumed. 
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APPENDIX A: Fortran Listing for Program LOBMUF 


The basic modifications to Blake's LOBEPLOT for program LOBMUF occur 
in Subroutine LOBES, as described in the Plotting Techniques section of this 
report. In addition, Subroutines PDSN, PD, INVERS, MARSUWR, DGAM, DEVAL, GAM, 
EVAL, and SUMLOG are incorporated directly from RGCALC. 
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STATEMENT 


PROGRAM L Sanur 
ReaALe@ TLA5L (10) 
COMMBN SYTPES 21,02,593,94,95,G6,ERROR 
C3MMONSADELTS/T ACD 
CemMMensCSC/ICSC 
COMIN SLT/LTS 
DIMENSTS™ 4(1500) 
PATA LT N/0/ 
C4LL RSSTCP 
READ 1, 4NTHGT 
ERRORSO 001 
C4uL PLATS (4,1599,0,5) 
PRINT 39 
39 FSR“AT (' Ta® DEN TURRET SBSITISNS 4RE USED, 010 AND D1}, 
wARTS wITH PEN TUPRET PASITIAN D190, CAARSER PEN [8 O11,'77) 
2 FEAD (5,5,8N03190) ILASL 
200 REAG 3, XMAX, VMAX, IMAXK yg HMAX,PHMIN, THMAX,AMFT,ROR 
NRORSRDR 
IE ¢ NROR LT, 1 ) GETS 1190 
SFac = Yvaxe, 125 
52 es ,175*SF ac 
Y & 1,5#SF4C 
TF (LIN EG, 18) GE TS ol 
IF (C¥¥AdeY),LE.9.5) GO TAH 61 
60 SFAC & 9,S/CYMAXeY) 
YMAX 3 VAN RSF AC 
XMAX 3 XMAXY#SFAC 
Y = YeSFAC 
os HeSFAC 
61 XsYe 3 eh 
CALL SFIGIN (X,Y) 
PRINT 10, ILSBL 
CALL RHSCHT (XMAX,YMAX,RMAX, MAX, ANTHGT) 
CALL LSBES (XMEX,YMOX,RMAX,HMAX WHET, THMIN, THMAXZNRDR) 
CALL LETTER (CO, oY pS, IL 48L,0,0,80) 
CALL GRIGIN (XMAKO3, #Y,eY) 
CALL REZER® 
GS T8 2 
100 CALL ENMPLT 
GSTS 130 
116 PRINT $20,%80R 


PLOT ST 


129 FSRMAT(//,!' PRAGRAM TESMINATIAN OLE TR ILLEGAL DATA ENTERED FAR NU 


*#MBER AF RADAPS , 
130 CANTINUE 
1 FORMAT (F10,0) 
3 FARMAT( 8F19,0) 
5 FerMar (1048) 
10 FRRMAT(3¥,1048A//) 
ENO 


o © « « NUMBER ENTERED [8 ',125) 


10 
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STATEMENT 


SUBRSUTINSE LARES (KMAX,YYAX,RMAX,HMAX, AHF T,TEMIN, THMAX,NROQ) 

THIS VERSI9 COMPLETED JUNE 1976 T* ALL4W PLATTING LA@RING PATTERN 
FAR MULTIFREGUENCY CASE WITH INCOKERENT INTEGRATIAN, 
COMMSN/CSC/ICSC 
CAMMSN SYTPES X22, Y2, X1L1e Yilso XA, VA, ERRSR 
COMMONSO/KEC1ALI,YYCL81) -CTLC181),S5N10161),0ELL 
CAOMMENSOUML/SPFF( 2000 ) 

DIMENSTAN FXNG(50) —NIM(50)4SN9(50) 
OIMENSTIN RFS(S9),FMH2(50),840(50),SL08(50) 
DATA TIMAX/2000/ 
Dara [csc / O/ 
Data PT/3,141592654/ 
DATA PT2/6,253185307/ 
Data RDON/,01745329252/ 
DATA C8NV/{ OUSTAAS3 IE CUS 
DATA 4€ / 2,786526684E7/ 
NROSSNROS 
RNSSFLAAT(NROS) 
RFAD 994,897, PFA,PULS,CASE,AMFY,TILT,PSL,CSC 
PRINT 905,°97,PFA,PULS,CASE,AMFT,TILT,PSL,CSC 
905 FORMAT (702%, 4mMPOTS,F§,2,39%,UnPF AR, F5,2,2X%,S5SHPULS3,F6.1,2X,5HCASES 
pF So, 2h, SHAME TS, FOL, 2X, SM TILTE,FS,2,2X%,4ePSLE,FU,1,2X%,UHCSCa,F 4, 
#1,/) 
NPULSSPULS 
KASE@CASE 
CALL POSN (PCT,PFA,NPULS,KASE,S0B) 
SNNZ10,*#(SO0B/10,) 
GAMMAS2,3/SNN 
RETAB1,/(1, GAMA) 
904 FARMAT (8F190,0) 
PRINT $90 
590 FORMAT(' THE CALLING PARAMETERS SENT TS SUBRSUTINE LSBES APF AS FS 
*LLSws') 
PRINT 591 , ¥MAX ,» YMAX , RMAY , HMAX » WHET , THMIN » THMAX, NFOR 
59, FIRMAT (AK, SHXMAXS,FU,1,2%,SMY4AXS,FU,1,2X,S5SHRMAXB,F6,1,2X,SHRMAXE, 
® FIL, 2X ,SHYHFTS, FS, 1, 2X, 6HTHMINE, FULL, 2X, GMTHMAXE,FUel p2X, SHNAC RS 
#,13) 
he @4 Trst,1IMax 
PFF(IIT) 2 0,0 
E Ss YMAXMRYAXS(XMAXSHMAX eCANY ) 
EXsxMAy/RMAX 
ZaWHF Ts, 3535534 
X¥P21,60e45 
YPa{,£eus 
Ce 801 JJ81,NROS 
READ (5,802,ENN3807) RANGE, FRES,BEAM,SOL98,PULNUM,FREF 
803 Icscscsc 
IPSL=PaL 
NUM (JJ) SPULNUM 
RFS (JJ) SRANGE 
FMH2Z (JJ) FREQ 


11 


arog 


redial 


STATEMENT 


) Aad (JJ) seb am 
SLOS(JJISSCLAB 
An2 FSRMAT (6F1U,)) 
SRINT 399 , \RDR 
499 FARYAT(/S/,20x,'THE RADAR PARAMETERS INPUT FAR RADAR NUMBER ',73 
a,’ ARE SS FYLLS4S') 
OEINT S92,5FSCSI) -AHETSEMMZ (SI) BMD (SI) »SLOBC JI), TILT, POL »CSC,PLL™ 
ely ,FREF 
592 FIRMAT(/,2X, UHRFSS,F6.1,2X,5HAKFTS,FO,1,2K,SHFMMZB,F7,1,2X,4HBnCs, 
F512 62%, 54SLORS,FS,1,2¢,SHTILTS,F5,1, 2X, UHFOLS,F2,0,2X,4HCSCB,F2,0 
ee 2X ,ORNPULSS,FS,0,2X,SHFREFS,F71/) 
IF CPREF NEO.) GOT? Sia 
CALL PASS (OCT,PFA,NUM(JJ),KSSE,S0R) 
SNRCIJ)B1G, ee (SOR/19,) 
. 68 T9 §15 
514 SAXR(JJ)SSN% 
RESC IJV SRFS (JJ) #SGRT(SORTCFYHZ(JJ)/FREF)) 
515 CONTIAVE 
PFSC IJ) SFFS (JJ) e(SQRT(SGRTCSNR(JJ)/§NN))) 
<RDORSNRDR ef 
any CANTINUE 
TILTR = TILT # RON 
TRMINRETHMIN@EON 
TRMAXPSTRMAX eRON 
CELLS ¢ THNAxX9 © THMINE ) f [Max 
PI2Z30E = 3,9*DEL1 + Pl2 
THETASTHYINReDFLI 
X 3 0 
INDEX = 9 
Fa1.5 
TOASrs1 
CALL CASHSAF. 
AR QSAHFTeAHFT 
HAE © 2, & OHET/(3,24E) 
AEH 3 AE w (SE © AFT) 
PARAM 3 SQHT (SE/(2,xAMFT)) 
98 810 ITs t,rI“ax 


PF4zi, 
TRET2 3 TRET? + DELI 
IF (THKET2 ,GE, PL230E ) GSTS 4o 


T2 3 TAN (THET2) 

$2 = SIN(TRET2) 

$3 3 $e 

PST s TRETO 

C8 ANS “s1,NROS 

RePsBuh (M) #PAn 

Bag 2 8aGe,5 

CSCT = ,7071 # SIN(Bw2eTILTR) 

TF (AnO(“),LE.45,.) GE TA 25} 
250 [REAMsny 

CANSTS9N,/E ad (™) 
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STATEMENT 


IF (Rem (4%) ,69,360,) I[S8FAYsel 

Go TS 282 

IBE AMS} 

CANSTS1.391S7/SIN (Awe) 

ANLL BARSTIN (PI/CSNST) 
CPATISI0,#4((13,20SL98("))*,05) 

* 8 GB3,S73/F“HZ(™) 

#4 39,25 * « 

wLIM 3 0,01 # & 

Fac = Plas 

PO 8 2, « AHFT # Se 

IF CINOEX ,E%, 1) GS fT" 77 

f2a3 = Tas 3, 

GAM 3 HAE/S(SCRT (T23ee2s4hE) e723) 
PSI = TRET2 + GAM 

$3 = SIN(PST) 

ZETA 2 PSRAM #@ T2 
DiSOS7735S#SURT(C1L ee, eZETA/SGET(ZETAeZETASS) ) 
POLSSORT (AH2¢AF RaGAMaGdM) ac, aS3eS3 
IF (ABS CPN@PL1) GE ~LI™) G3 TS 79 
IF (C1,L7.0.999) G# TS 79 

01 31, 

INDEX = 1 

FD = PD} 

IF CPO LT *4, 88, THETA. 7.0,000873) GS TS 300 
G8 TS& 30% 

IF (ICASH,NE,1) G8 TS 300 

CALL CaSRSF 

ICASK30 

CANTINUE 

CALL SEASFFE CFMHZ(M),PST,TPERL,RHS, PRT) 
PIZSPTeZeS$3/a 

PIZ2sPTZeP IZ 

RUF SEXP (oA, #PTZ2) 

ANG 3 THET2 @ TILTS 

IFCIBEAM) 269,261,262 

PaTsi, 

Ge Te 24 

PATSCOS(CENST#ANG) 

Ge TS 24 

IF CICSC,NE.1) GO TS 61 

IF CANG.LE.8*2) G* TS 61 

Par = CSCT/S2 

Ge TS Q4 

IF (ANG,NE.9,) GE TS 23 

Par 3 {, 

Ge Te 24 

UU ® CANST # SIN (ANG) 

INT & WUsPT2 

CIFF 2 Uy © INT # PI2 

TF (ABS(4NG),LE,4NLL) GE TS Lol 
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STATEMENT 


CPAT 3 CPAT] 

GA Te {62 

cPart 3 j, 

Part 3 CPS4TRSIN(DIFF)/uUU 

4NGR 3 TRET2 + TILTR # 2, * GAM 
IF CIBEAM) 270,271,272 

PaT@si, 

GS3 TS 27 

PATRECSS(CINST#ANGR) 

Ga TS 27 

TF (CANGR,NE,O.) GS TT? 26 

PaTR = 1, 

G0 Te 2? 

UUR = CONST # SIN(ANGR) 

INTR & UWUR y PT2 

DIFFR = UU® ew» INTR # PI2 

TF CABSCANGR),LE SLL) G9 TS 66 
cPpart = CPAT{ 

G"® TS 67 

cPar sl, 

PaTR = CPAT#SIN(DIFFR)/UUR 

IF (ABS(P4T),GE,1,Ee44) GA TS 29 
IF (P4T,GE,0,) GS TS 281 

C 3 eRHO RUF ePATR ae] .EUSeD{ 

30 TS 30 

C01 eRHSeRUF ePATRe!] EUS 

68 TS 3¢ 

D s Diw PHY w PATR/PAT «2 RUF 
ALPRA = FAC «# PD + PHI 

INTL & ALOMAsSPT2 

CIFF{ £ &LPHae [NTI ePT2 

Fo s ABS(PAT)& SGURT(1,.0 + DeD + 2,02DeCSS(OIFFI)) 
RXNG(M)S(RFS(M) aF ) wed 
PFUSPFUSSIRT(SGRTI(RXNU(M)))/RNA 
CONTINUE 

RNGSPF4 

RNG2ZER'GHRNG 

RNGUSRNG2 RNG? 

FACMINSIO, 

SUMEN39, 

Df 996 “Ms1,NRDS 
RRATIOSERXNG (YM) /RNGU 
PNOSFLAATCHUM(MM)) 
SUMENS(PNAS/PULS) #(RRATIO/( 1, +GAYMA/RRATIC) ) ¢SUMEN 
CONTINUE 

FACMN2sF AC™IA 

FACMINSSUMENCBETA 
FSUMSEACYNZSFACYIN 

TF (CPNG,ES,0.) GS TS 993 

IF (48S CFACMIN),LT.0,91) GE TS 993 
IF (FACYIN,GE.et.) GE Te 902 


ee ieee nS sh alt 8 


Por ed 


90s 


902 


993 
999 
aya 


ent 


2 


STATEMENT 


maNU4BSATS TOF ACMIN) 
PAC’ IN SEACH TN eHINUY 
CONTINUE 

RPREVSRNG 
RNGGSRNGSe (lL oF dCIN) 
RBNGSSGET(SIRT(SNG4)) 
ROIEFZARS (ROREVeRNG) 


TE CPOTFF LT, O.t.4NC.FSUY.LT.0,91) GA TE 993 


6° 78 997 

CANTIAVE 

PFFCIITJSRNG 

CANTINGE 
TRET2ZETRY[ NR eDELL 

fe 896 TIS1,7I“4x 

AsPFP CII) #E& 

TRET2 = THET2 + DELI 
NaNet 

Xs 4 a CIS( THETA) 

y 3 Aw & #& SINC THETA) 
TE CX LEXA) GS TH SO1 
TSu“¥88 

G* Td Sv2 

TSUMXsn 

IF CY,LE.Y"4u) GS TR Sha 
ISsu“Ysa 

Ge 9 505 

TSU¥YS0 

IF (YP LE, XMAX) GA TS SOT 
ISUMXPs2 

6s TS 5998 

TSUMXP30 

IF (YP LE,yMax) GA TS 510 
TSUMYSst 

Ge T8 $11 


9 TSUMYSs) 


TSUMBTSUMXF+ISUMYSTSUMXP eI SUMYP Ss] 


BE TS (6914402, 60 5,602,605 7620 4629, 620,609,620 ,620,620,605,629), 


*Isu¥ 
IF (N,NE.1) GOT 3 
CALL PLAT (x,Y,3) 
x2au 
y2syv 
yPax 
yPsy 
Gere @n6 
IF (N,NE,2) GE TS 440 
yas 
yasy 
1Psx 
yPsyv 
xLi1sx 


15 


t 
‘ 
a 


400 
620 


692 


621 


605 


009 


an 
650 


651 


652 


STATEMENT 


YLLsy 

GST" A. 

COLL MINTAP (X,Y) 

xPax 

YRsY 

IFCTRET2,GE,THMAXR) GI TS 40 
GSTS Bix 

TF (N,ES,1) GS TE 620 

CALL OSHS 

TDASHZ0 

IF CIAXIS,EG,1) CALL PLATCXMAX,YMAX,2) 


COLL INTRST (Oy, YAK, eMAK,YYAK,X,Y,¥P, YP, 19, YC) 


CALL PLATECXKO,°Y0-2) 

CALL CASHOF 

X2sxo 

yasyo 

XAsX 

X113X 

XPsx 

yasy 

y1isy 

YPsY 

Ns2 

GST® 806 

TE (N,EG.1) GE TS ean 

CALL DASHAN 

TDA4Sks0 

CALL JNTRST (XMAK,0,,XMAK, YMA, X,Y, KP, YP, X0,10) 
GS TS 621 

IF (N,E9,1) G9 TS 620 

CALL MINTAP (XP,YP) 

CALL PLATCKP,YP,2) 

CALL INTEST (0, ,YMAX,XYAXK, YAK, K,V,XP,VP,X0,Y9) 
TAxtSs2 

CALL PLST(K0,Y0,2) 

Ge 73 629 

TF (N,E9,1) GO TS 620 

CALL MINTAP (xXP,YP) 

CALL PLOT(XP,Y9,2) 

CALL INTRST CXMAX, 0g XMAX,YMAK, Xe Ve XPLYP,X0,Y0) 
TAXTSst 

CALL PLOET(KO,Y0,2) 

Ge TS 620 

TF CISUM,NE,1) GP TP 651 

CALL PLET(X,Y,2) 

RETURN 

TF CX ,GT,XMAX AND, XP,GT, xMAX) GA T*® 652 

Ge TS 653 

CALL INTREST (XMAK Oy XMAK YMA go Dg pKa V eK V1) 
CALL DasHe™ 

IF (YO,LE,Y'Ax) GO TS 655 
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“ge sh matali Rt a 


STATEME''T 


Ss CALL FLAT (rsx, van, e) 
CaL& Purkgr CO VOM RAK YMA Cg De Xe Vp XO V0) 
o55 CALL PLST(¥I,¥C,2) 
COLL CF aSrse 
RETURN 
SSS IF CY GT yuan AND, YO GT yMAx) GF T® 655 
Ge T* 657 
656 CALL CaS4A® 
CALL INTHST (Ope YMA, KNAK,YMAX, 0g p0 ee Xo Ve X00) 
IF (xO, LE,%“4x) GF TS 059 
S27 COLL INTRST CxXMAK OL XMAL,YMAN Dg eS g yp Kap X0,Y0) 
boo CALL PL*T(x0,Y9,2) 
COLL Dass 
aS7 FETURN 
659 IFCIAxITS.ES,1) CALL PLATCKMAK,YMAxX,2) 
6e T® 666 
BGO CANTINYE 
GSTs ay? 
ROT PRINT ace 
ROR FARMAT({Ht,! END PF FILE BAS BREEN PE4D “RERE CATA CARD SHOULD SE!) 
ang ENO 
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Spee ee _~ 


ee ae ae ee SOLOS rye oe Tee eae 


~ere STATEMENT 


SLBRSAUTISE BOS. (POT, PFA, FULS,KASE,S98) 
J} FYTERMNAL pe 

DIMERSTIS NB0C5)-SLSPE (CS), POR AC (5) 
Cavvsn sSEOS/ PFLAST,NLAST,KSLAST 

DATA PELAST/N,/ 

CATA NLAST/O/ 

DATA K§LAST/e1/ 

CATA PRUs 2 Set deol dent Se2vtSee/ 

DATA SLSPE JO, pT epBee ter ha/ 

BATA POF AC/E B20 gp QV eel beni de/ 

DATA CEYIN/@30,/ 

DATA OSM4X/50,/ 

MATa CRLS9,S 

DATA CS2/0,/ 

Data POLAST/C,/ 

TF (POT. NE.POLAST) G* Th 2n 

IF (PFA,"NE,PFLAST) G& TR aU 

IF (NPULS,%E.NLAST) GO T* 20 

IF (KASE,ES.*SLAST) RETURN 

PRLASTSsPCT 

PFLASTSPEA 

SLASTENPULS 

KSLASTSKASE 

K3K ASE + { 

PULSSNELILS 

PRLSOBO (HK) OSLIPE(K) @ALSGLO(CAULS) © (PDT|,5) ePODF AC (K)o(PF 408, )e uel, 
DRe2s0R1 +2, 

C4LL It VERS(RAMI®S 7CRAYAX,F81,°82,4,15,N81,S8CR,P901,FP5T, Pr) 
ENS 


2Fintwe 


cope tear, aramrreceisinimnsneistas ee . a 
hee ™ a a a es eee aera ee Menem ee . 


STATEMENT 


FUNCTION PFOSKDOS) 
CeMYSRSOOG/SEA,* p< OSE 
‘PsN 
FANSFA 
KASSKASE 
° CALL MERSAS (SNDR,NE, FAN, KAS,POL) 
FPOSFDI 
EAD 


} mit DY 
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ee Ce CO. Oe. | 
be ; pactetiteneeoen inet ine Seni 9 Te Andee 8 meet 86 a . ; 


STATEMENT 


SUBRPLTINE INVEFSCEMIN, MAX KLM MAT,NSIG LIM NOT KF LS FTL) 
TEST = 16, ee(erSlG) 
Fr 3 FT 
If CET EG, 0.) FO @ ft, 
NOT 2 1 
DELTA ws ae] © x4 
XY sxe 
X23XeT] 
Fisk (x4) 
F2sF (x2) 
SLOPES(F20eF1)/(x20x1) 
TE CSLAPE,NE,0,) GO TA 21 
to FMAxsR(x4x) 
PMINSECXMIN) 
SLOPE = CFMAX@FMIN) S(XMAXKOXMINY 
2! IF CCE 2eFT) eSLOPE.GE,9.) GS TS 23 
22 xt{sxe 
Fisk2 
K2sxX2e0FL TA 
TE (x2 .GT, xXMAxK) YoRKvax 
FasF (x2) 
Ge Tf eal 
23 TF CCF TOFI1) *SLSPEGE.0,) GI TA 25 
Pu Kean] 
i KisxLeMELTA 
t TF CMY aL T, XMIN) XPSHMIN 
Fask\ 
Fisk (x1) 
69 TS 23 
25 xAsxt 
xBsx2 
j FQSFQ 
; TF (OR§ (F2e0eFT),GE,ABS (FISFT)) GE Th 6 
7 F2esF2 
FasFi 
FisF22 
X223X2 
Kesxt 
X}sx22 
G8 T9 6 
1 Fysr Cx) 
x{2x 
TESTI s ABS(CFISFTI/FO) 
; TF CTESTIL,GT,TEST) G* T® 6 
| 2 RETURN 
& [PONCE LY.GI4) GO Te Us 
12 PRINT 40 
PRINT ul, LIM 
PRINT ud, XMIN, XMAX, XL, XHU, MSIG, LIM, NAT, Ky Fi, FT 
: RETURN 
i 13 IF CF1,NVE,F2) GS Tt to 
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ee on er 


“ew eal Hh, e 


STATEMENT 


f 
1S IF CEL vE,e8) Ge Ts ya 


17 «Ssxt 
xy251F 
89 T& 49 

/ {= wasxt 
F x2sud 
ij 19 xa(thoxd) 2,5 
E G9 T9 3 
i 1> XS(KLex2)e(FTeF2I/(Flers) + x2 
E IF ¢¥ LT, x4) vod 
i IF (4 ,GT, “«8) s25K8 
a NSE s Nay © { 
B FosF) 
. 122x1 
c 60 & 4 
a FSEYAT (7s! MESSSGF FROM™ SUBERUTINE INVEPS ee '//) 
& wo} FRevar (1 FUNCTIAN IN vEGSTSt NAT ACCOMPLISHED wITMIN SPECIFIED ! 
t _* o15, ' ITEGATIANS,'/) 
Ke 42 FSRMaT (ys! INVESS FATA ETERS WERE ',4(E10.3,2),3(13,2X),E10.3, 
B . 2fex,Fy%,$) 4/) 
) AC 


| 
kK 
| 
| 


~ 
a 
é 
2 
a 
. 
tz 
3 


E 
r 
' 
| 


4 
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or 


3 
@ 


uy 


ie hae a a ee ee ae a ee meaner Per er ee. ae . — . i i, “e 


STATEMENT 


SURRfUTINE MBPSwR (SNOR,%,FA4,KASE,PR) 

PeUBLE POFCTSISS EXNPE,YHPP,GAMER,PYR,F,YU,E0,YL,EL,STEP, YS 
DAURLE PRECISISN NGO, DEVAL, SUMLOPS, SUML, FAN, F* 

SAR & 19, e%e(SNDA*e,1) 

“ODES 

IF (MEME) 400, NO, FOO 
FANSOLSGLQCOLIGCUSCOISOLAG (Lo (10,00) ea(oFA))) 

GS TO 995 

FAN s FA 

IF(N) 99,99,2 

IF (F4)99,99,3 

IF (K ASE) 99,4,4 

IF (KASE ed) 5,5,99 

ENPS = Qe 

ENPR s FAN 

EN 2 N 

YBFR = 9, 

IF (NPREV ,E€G, ' ,AND, FAPREV 
IF(Nw{2) 7,7,8 

YRAPRSEN&( 1, s2e2eENPR/EN HH ((2,0/3,00) 40.01 5#ENPR)) 

GO Te 11 

YBPR S&S ENa( LoL yg SeEPR/EM we (54,011 ENPR) ) 

ENPR 3 10,%8FNPR 

GAMPR s DGAM(YSEPR,Ne1) 

PYR S ,See(1,/ENPR) 

SUML & SUMLSGCNel) 

ITF (GAMPRePYA) 10,12,12 

es att 

G8 Te 14 

mos et 

YRPR 

CEV4L (vO,Ne},S'"L) 

VooH ! 
CEvVAL(YisNel,SuML) ; 
STEP = GAMOR +» He(ECSE1L)/2, | 
TE COCSTGN( 1 DO sSTEDSPYROOSTON CL OD H I gE Te0e) GR TT 20 
Yo = Y1{ | 
Eo = Ft 
GAMPR s STEP 

66 TS 14 

IF(*) 22,246,246 

YE & Y1 © Me(PYSeSTEP)/(GAMPReSTER) 
Ge T3 30 

YR & YQ + HH (PYReGAMPP JS (STEP eGAYPR) 
83aS 3 Ye 

YR = ArTAS 

SPREV = ': 

FAPGEV = Fh 

y s §'e 

K 3 KASE4] 

GR TS (106,290,300,409,500)6 


,&0, FA) GO 8 777 


< 
- 
oe ot ae at 
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STATEMENT 


me SUM za 
P s EXNex 
TF (vRer eet) 150.1972,192 
eR vets )/2q © CSISTECCEN@{,)/2,) ee2eP av) 
“and (4Se%) 
1 eGAN CYB, KSeNet, Th) 
Eva_(P,«kS) #6$ 


x 
o 
oo 6b 68 oof 


TE wT 
110 TEMP = SUMSTERY 
TE(CSUMeTEMP) 112ptloette 
tt2 Suh = TErP 
TECK) {15e116-114 
tia TERM s TERMePL IST (K)e(GeTLI/(FaG) 
G = Get 
= Ae] 
2 TLePLIAT (KON) /YE 
te Mo 


1146 TN®YRSFLAAT (KS4N) 


TEQM & TSePeG/(GSeFL 74? (k)) 
Ee 129 TEYP = SUv+eTFRé 
. IF (CSUM%eTEME) $22,190,190 
122 Su” = TEP 
TL S TLeYS/FLGAT (Ker) 


« = key 
- TERM = TER “wee (GeTL)/(GeFLIAT (K)) 
G§ 3 GeTL 
+ G9 T° 120 
4 150 KS 3 of, © EN/2, + OSGPTCENSH2/4,4P HB) 


KS 3 "AxkY (KS,0) 
GS 2 Gat(YA,kSevej,T) 
IF (GS) 174,174,155 
155 TS = EVAL(P,KS)*#SS 
G6 = GS 
TERM = TS 
kK s&s 4§ 
TL © TN 
160 TEMP = SUMeTERM 
F IF (SUMeTE4P) 152,166,166 
. 162 SUM & TEP 
IFCK) 165,166,164 
164 TESM & TERYRELOAT (K)#(G+TLI“ (PRG) 


! & 2 GeTL 

. EE = TLeFL Sar (KeNow{)sy8 
i € § «ef 

r 69 TS 192 

‘Ss tho TL 2 TNRYR/FLOAT (KSsN) 
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Sa Pape are 


Pa or 


RE Ln ee Se ee e.g ngs se ah mont AR eS 


STATEMENT 


K s KSel 
G = G$eTL 
TERM = TSePeG/(GSeFL IAT (K)) 
170 TEMP s Slit # TERM 
TF (SUMeTEMP) 172,174,174 
172 Su” s TEMP 
TL 2 TLaeVAsFL AAT (K+) 
TERM s TERMePae(GeTL)/(GaFLIAT (K41)) 
G = GeTL 
K Ss Ke} 
Ge re 170 . 
174 SUM 8 1,eSuU™ 
190 PN 3 §UM 
6¢ 78 99 
290 ITF(Net) 210,210,220 
210 PN 8 DEXP(eYB/(1,+x)) 
ge To 90 
220 TEMP 2 1, © LesCEN#X) 
; PN Ss 1, © GAM(YB,Ne2,0U") © DEXPC(ENe1, PAL AG( TEMP) eVB/ (1, EN eX) ) 
* 2GAM(YB/TEMP,N@2,0UM) 
Ge re 90 
BOG TF(Net) 310,310+320 
a $10 PN &S DEXP(eY¥B/(1,%X)) 
; G® T8 99 
320 PN 3 4, © GAMCYB/(1,4X),Nol, DUM) 
Ge Te 94 
400 IF(N@2) 410,420,430 
“ 410 PN 3 (1 ee2yeXwYB/(X42,) e022) HNEXP (2, #¥98/(2,4X)) 
6e TS 90 
420 PN & (1, eYB/(1,%X)) #DE XP (eVYB/(],4X)) 
Go T*® 99 
Cs 2,762, +ENex) 
0 8 1,°C 
IF (YBe#DeEN) 440,450,450 
auc SUM = 0, 
TER“ si, 
J @ N 
4u2 TEMP s&s SIM+TERM 
IFCSUMeTEMP) GUU, 446,446 
uuu SUM = TEMP 
TERM 3 TERMweYBeO/FLAAT (J) 
i J 3 Jel 
GO TS uue 
i 446 PN & 1, © GAM(YB,N@2,9UY) © CeYSwEVAL (YB,Ne2) 
“ * + DwEVALCYB,Netjyw(l,eCeVRe(ENeo2,) #C/0) «SUM 
; G6 Te 90 
i $50 PN & 1, © GAM(YS,Ne3,0UM) # YOREVAL(YS,NoS) eC /C 
7 * + DEXP (eC eYSe(ENod,) eALSG(C))e( le sCeVBa(ENo2d, )#C/0) 
q. * GAM (YBel,Ne3,0UM) 
ws Ga TS 90 
vs 5co SUM 8 4, 


ae ow 
c 
uw 
2 
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SS Se a 


a — —_ arn - alk b kke wm alae Neh - ; Sax es ses 


% 


a Pp eat joe 


A I i iene: frp cent in i i a ail 8 


STATEMENT 


2,/(2,**) 

1,eC 

Csh 

Cev3 

KS 5S C3, HEN O(YEUD))/2,e0SGITCCEN@L He (Y90D)) we2s/d,o(YReD) #(ENOl1,)) 
K§ = “SIMO (h#§,*) 

K§ Ss “AYO (KS,%) 


Veou™"™!] 
oo a os oe 


FK§S 2 kS 
kK S MINQ (KS,N) 
TEC YBer Ne (1,49)) 550,891,501 
Set GS S 1, = GAM(P, 2 el ox§,T*) 
IF(GS) 526,526,592 
S02 TS = DEXF CF KS AL SG(C)+CEN@FKS) #4L96(0) +SUMLAG(N) @SUMLOG (KS) 
* @SUMLOG(J)+4L96(G6S8)) 
G = GS 
TERM = TS 
TL @ Th 
Sic TEMP 3 SUM*+TER" 
IF (SUMeTEMP) 512,516,516 
512 SUM = TEMP 
IF(K) S16,516/514 
5146 TL & TLeP/FLAAT (2ehek) 
TERM & TERMHELSAT (K)w(GHTL)/(GWPLIAT (NeKo1) 2G) 
G s GeTL 
« S Kel 
6@ TS 510 
S16 IF(KSe") 518,526,524 
S18 TERM & TSetaFL@AT (NekSpae(GSelV)/S(CFLSAT (KS01)*#68) 
Gs GSef* 
TL 3 TNFLIAT (2#Nof[ekKS$)/P 
K Ss KSel 
520 TEMP = SUM*+TEQM 
IF (SUMeTEMP) 522,526,526 
522 Su“ s TEP 
IF (KeN) 524,526,526 
S2u TESM & TERMaGeFLIAT (Nek) we(GeTLI/(CFLIAT (K41) 2G) 
G s GeTL 
TL S TLeFL AAT (2eheoejfok)sP 
K Ss Ke] 
68 79 529 
526 PN & SU” 
68e T8 99 
550 GS s&s GAY(P,2eNeo{ek§,TN) 
IF (GS) 576,576,552 
$52 TS Ss DEXP(FKSwALSGIC) O(ENSFKS) eAL SG (0) *SUYMLIG(N) @SUMLIOG(KS) 


* eSUML*G (J) o4L°6 (68) ) 
G 3 GS§ 
TES” 2 TS 
TL & Ts 
25 


~~ 


me 


Son 


Sa2 


NPBA ln ie arin tant ye a 8 sk mln A 


STATEMENT 


TENE = ZlveTERY 

IFCSUM“elENE) S6e,S660566 

SUM = TENE 

IF (KR) 36m-5667564 

TL S TL POEL GAT (2e\en) 

TEAM = TESMePLOAT (K)e(GeTLI/S(CePLAAT (NeKe{ ) 9G) 
G s GeTL 

kK 3 Ke{ 

Be Te San 

TF (KSe%) 569,576,575 

TERM S TSeGeFL OAT (NeKS) a(GSeTRISCPLIAT (KSo1) 26S) 
& = GS§efi. 

TL f TNRFLSAT (Atel ekKS)/2P 

« Ss «Sel 

TEMP = SUM*eTERM 

IF (SU“eTE“MP) 572,574,576 

SUM = TEKe 

IFC KeN) 574,576,576 

TERM s TERM eM ePL AAT (NeK)a(GSTLI/S(CFLSAT (Ko1) 4G) 
G = GeTL 

TL = TLeFLSAT (2ehe{ox)s/P 

K s Ke} 

68 TS $70 

PN f {,0Su™ 

65 T93 9) 

IFCPN) 91,94,92 

PN Ss 0, 

G3 Te 94 

IF (FNw1,) 94594,93 

PR 3 1, 

DETURN 

WRITE (81-9) “FA ,SN&,KASE 

FSPMAT (140 /504 UNREASANARLE CALL SEQUENCE T& “ARCUM, 


ZERe RESULT 


74S GIVEN //4H N & T8,Sx,SHFA = E1648,SX,5HSNR & 
E15,.8,Sx,6rKASE 3 ITA) 
PN SQ, 
RTagS s ), 
RETURK 
END 
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~Se 


STATEMENT 


FUNCTIAN OGar¢(e,") 
CeURLE PFECISTS:: Siw, TEM, TEMP,FJ,0GAM, SEvOl, Fy, StL, SuYLAG 
SU” 3 OG. 
K se 
TF (KON) $0,200,209 
oo J 3 NoY 
j SUML = SuUMLeG(S) 
fe TERM = DEVAL (S,J,SU"L) 
: 16 TEMP @& SiMeTERM 
bt TF (SUMeTEMP) 15,29,25 
. {S$ SUM # TEMP 
4 Jz J+1 
4 FJ 3s J 
TERM & TERMNHRSES 
Ga T® 10 
20 CGAY = SLM” 
q RETURN 
g 200 J se N 
SuML = SuMLIG(J) 
t TEP“ & DEVALC8,J,SUML) 
% $0 TEMP & SUMSTERM 
TF CSUMeTEMP) 35,40,40 


et 35 Su” 3 TEMS 
o TF (Je1) &6-36036 
is Jo Flag 
; TERM & TERMwFJs& 
# Ja Jet 
. Ge Te 30 
gl 4o DGAM = 1,°Su™ 
if RETURN 
4 ENO 


ye 


PS ee RAE 
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er ee eee ee ~~ 
ee eR 


STATEMENT 


FuMcTyTes DEVAL C¥,%e/SuL) 

OMURLE FRECTSTIA®® cPAXN,EN,CEVAL, Y,SUML 
xOQ" = @Y 

IFQN) 29,20010 

EN 2 N 

KPSN & KPANSENSDLIGCY) @SUML 

DEVAL = CEXP(XKPAN) 

RETUGA 

ENO 


wer % 


vm ye sRaP 


ae 


eet ee 


STATEMENT 


FISCTIA® GAM(A,N,TA) 
SLY 3 0, 
« 26 


TE (Ken) 1UU2 2052209 
Jos “of 

TEPM s EvAL CR.) 

TA 3 TER“eF LAAT CJ)SF 
TEMP s&s Site TER™M 

IF CSuMeTE4P) 15620,28 
Su» = TE® 

Jo= Je! 

Fl zJ 

TERM s TER’ RRSP) 

Ge Te 44 

GA“ = Sl” 

RETURN 

J sk. 

TERM s EvAL(8,J) 

™ = TERM 

TEMP = SUMeTERM 
TFCSUY%eTE"P) 35,409,450 
SuY & TEMP 

IF(Je1) 40,5056 

FJ @ J 

TERM = TERMePJ/B 

Js Jel 

6a Ts Bo 

GAM Ss 4,<«SU™ 

RETURN 

ENO 


Le AA EAMG he RN AB AGE 9 sg PL 


—a —————— ideal \apeappen alc. ik cael li a 8 i TH 
any 2, - 


STATEMENT 


FUCTIAN EVAL (Y,') 

XPEN 3 eY 

IFON) 20,2°.10 
19 FX @N 

KPO S KOSN SEN HALPG(Y) @SUMLAG(A) 
2¢ EVAL 3 ExP (xPRn) 

RETUR?. 

ENO 


OO Pa IEE 


STaTEre st 


PL E?LTan Si 4L*G6(%) 
MSPSLE CPECTSIS. &, Fe, SUMLAG 
NIVMENST®* SC1000) 

ROTA Cuts £ Oes 

Carvh CumKe ¢ OF 

Parva AL4astsys 

AM ARSE GO 

IF ¢Gusae vs) 2leldsed 
ruwa st, 

Cums g ht, 

uLast s 1 

4(1) 3 4, 

we 3 JOSS (h) 

TE (NSN ot) gu, S00 SS 
Su¥LSG 2 9, 

PETUGN 

TF Ch NeMLAST) 50,550,009 
SUMLSG & AMHR) 

RETURN 

K 2 LAST#L 

TE (NN@NMAK) 79,759,480 

Ce 72 yan, 

4(I) = A¢TeL) + DLOGCOFL*aT(I)) 
SCAST @ NN 

6s 8 §9 

IF (NL AST eN4AY) &9,90,90 
Qe Ay Tau, aXx 

ACI) = A(Tel) + OL®GCCFLYAT(I)) 
NLAST |= AX 

Bos A(x) 

kK Ss NMA Ot 

C9 92 [set 

Ros 6 + OLSG(OFLS47(1)) 
SUMLOG = 8 

RETURN 

ENO 


31 


| 
| 
| 


s90 


ast 


492 


agy 


u94 


u95 


51 


STATEMENT 


SUPRAUTINE RRACHT(KMAX, VMAX, RMAX, HY AY, ANTHGT) 
EaTERN aL Ft 

PIMENST® SNCLAL), R°GICLEL), TNICITE1), Jaleo) 
RIMENSTON K27(61660),¥26189) 

fIMENSTS* TREAC (5), IHF AC (4) 

CAMMSN sds REF, GRAD, FOF, C*X“ST, UC1E2), * 
CAMMON SMTPES K2Q2,Y22e01L Vile XAA, VOD, ERRAG 
DIMENSISS T4'G(9),4a'GC9) 

COMMON SBS XXCL4AL)» YVCLSH), CTLCLS81), Shrcretd, DEL 
QaTes EFRSR / ZOULS 

MaTa REFS, 20031 35~/ 

PATA GEANS,QONNK3RLR AB, 

DATA TANGA, 20,31 pSi ep Tt tare t5i ei 7e,t8i/ 
DATA BANG/D Gp Soler derSerlar dare p Me/ 
PATA THEACS§5 C190 SOM LN0O SHI L IODC 50006 LCIMOS 
DATA [AFAC/5,10050,19925007 

DATA JO/S3L Stele eteiSios%as 

CALL PENCHG(19) 

MELSXMAKe CY 

IF (CYMAK,LT,XMAX) CELBY"4X #04 

BA S BOTH LISS MRMALHEYYAXS(MMAX HXMAX) 

BA2 B SAvhs 

PREC & ,f Onn, 

CSNST = ,3048/1852,9 

RPADS2HEFRISO  LRSANTHGT 

AB 3 1,0 SEF 

AB2 = ASe4é 

CO 3 2,0 * REF + FEF wREF 

ELEV = e,ve2 Pi 

Ifs5 

NA 29 JLS1,6 

GA TS (390,491,092,4593,494,u95), JI 

40FL = ,%2 

MMs26 

Ge Tf? §1 

S0EL = gt 

"msgs 

Ge Té §1 

4DELS,5 

“ms20 

63 T8 §1 

ADEL = 1,90 

“325 

Gs TS Ss} 

ADEL32,5 

““ei2 

Ge Te §1 

4DEL 35,9 

Musy 

fe 29 In = 1,%™ 

TIslie} 
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NS i a te = = ANP Mw scat . slid — a ae ees one. sch ete tha th Pe ee 


: 


Seu 


Se} 
29 


29") 


STATEMENT 


ELEVSELEVS4S0EL 

BQ S ELE V/S7, 29577957 
SNCIT) & SINCROM) 

S$ s SN(TI)*#2 

LCIT) =s 432*5eC9D 

TF CII.SE.18!) GS TS Set 
TN s TAN (RON) 


TNY CII) = Tar, (RO) 
CTLCIT) = CASCFONL) 

SNICIIY = SINCRONL) 

Ge Te 29 


ETHCIEY 2 0,0 

Suictl) = 1.9 

CONTINUE 

HeRAwxax 

C3 e9¢ [xs 1,180 
RATSSGPT(CeA2e(TV1 CI) ) #82) 
XZCIM) Sava er asesT 
YZCIX) set i C1x) /B4T 

“1 3 0, 

Ty4% # 109 


IMax = INT (ALCGLUCHYAX) @ ,477122) 


MINTSIC, ##(J4AXe]) 

TJD & (MALS (1G ee (SA kot) je 
bexslJGrlder(JMAve)) 

68 30 J S 1, JAX 

IF (J be. J“4x) Ih = ISG 
TF CJ ghEst) GY VM F 

iy os £9 

69 T®8 a 

TJ = es 

De 30 1 = IJeT 4,19 

we wT # 10 ew (J = 1) 

TF (M2,EG,H4X) CALL PENCKEC11) 
C@ 304 K S&S 1, 38] 

Ns1é2en 


0991) 


FF CH2 sGdy 19 5) BsG6i(%) 2 0,0 


TF (HAE 3,19 AND, E81) GO TS 6 


GS t€ 7 

GAM s REEeGwacsag 
RI 8 1,0/RAD 

GG 2 2;0#(RI © GAM) 


GNG2 & (CINSTRAB/GR) #2, 0eS09T(CGeH2) 


Go Te A 


CALL SINCS% CHE,Me,PFEC ,15e8TC,hO1,%,F1) 


RNG2 © 2G1C) * RINC 
RNGI(N) & FNG? 

IF (he LT, HINT) GY TR 304 
A Ss RNG@exV¥ax/Orax 

R 3 Bhed 


33 


STATEMENT 


TF (NOME WSL}? GO FO a? 
66 CALL PLOT (0,,8,3) 
Xx(181) = 0, 
Yy(161) 3 2 
XS0, 
Ysa 
x22s0, 
¥2238 
G6 Te Zug 
a7 sLAST = xX 
YL4S? 2 Y 
MobaexZ(%)/KMAX 
YsSe¥Z( +) /X" ax 
aa TR (h2,NE,HMAX) GO TS 68 
&7 xX (N)s¥ 
YY(N)sY 
68 IF(K,NE.2) GS TS 74} 
780 xiisx 
Yiisy 
xaaax 
yas 32 Y 
69 79 304 
7A, TF CX LES CXZIN) %,0001)) GA T% 763 
THS CALL INTSEST CXLASTL YL AST XeVOX7ONIAVZ(NS AZ IN OLD ,YZ(%O1),x0,%0) 
COLL MINTAP (x0,70) 
CALL PIST (X9,Y),2) 
TF (H2,E9,r“%4x) GO TE 7AaS 
G8 Te 395 
7&3 IF (K,ANE,1481) GP 79 7e8 
i 7A7 COLL MINTAP (X,Y) 
f C4LL PLAT (xX,Y,2) 
XCSR s « 
yYeeR 3s ¥ 
G8 TS 300 
786 CALL MINTAP (x,Y) 
304 CONTINUE 
305 iL & 62 
0 CSS TINUE 
Ge Te 789 
‘ 78S CALL PLIT(XZ01),0.,3) 
CS 786 Nest, 
MM CNK)BKZ('K) 
YY(NK)SYZ(NK) 
TRH CALL PLOTCXK(NK) YY (1K), 2) 
CALL PLIT(K9I,YO,2) 


xC@R s&s 40 
YCSR s Yo 
f 789 x30, 
é Ys yMax 
- CALL PLIT(X,Y,3) 
* Ys 0,9 


34 


STATEMENT 


CALL PL IT (CK, Ye 2) 
aSx+ex 
CALL FPLITCK.¥,2) 
E COLL PENCHE (19) 
KAMERMAXOL, 
INT&S1A 
k TFCRMAx LT, 109.) INTRSS 
i IF CRMOX,GT, 3O0,) [in TREa5 
C8 31 KASINTR,KAY,INTR 
&G = KA 
Bs RGaex’axsPrax 
j BR s SAed 
Es Ne 32 «OC Ss te fet 
eosKkC 
IF (KC,NE.1481) GS Te G2 
91 CALL MINTAP (0,8) 
; CSLL PLATCO.,t,2) 
i Gs TS 31 
G2 KsdexZT (KC) Seman 
ee TF (X,GT,xx¥(1)) GE TT? 31 
y vshey7Z(<C)/xvax 
° TF CY,LE,CYY (KC) %,2091)) GA TH 72 
3 73 TF (KR,NE.*.O1) GE TS 731 
a 730 x} sxC9r 
YIsvCr 
Gea Ts 7$2 
a 731 xpexx(e Rot) 
‘ Y! = YY(KReL) 
; 732 K2SxK (KR) 
& Y2 3 YY(KR) 
é XAsx 
# Yes 
4) CALL INTRST (X1,V1,x%2,Y2,X4,V4,xXR, VE, x0,V0) 
2 CALL PLYT(KO,Y0,2) 
£6 T3 31 
72 vk sx 
yas y¥ 
IF (KC,NE,1) GE TS 86 
&5 CauL PL°r(x,y,3) 
x2e=nx 
Y22sY 
t 63 T9 38 
a : be TF (KC,NE,2) GS TS Bol 
“ Boy ¥L{se 
yijsy 
KAASX 
j YAAsy 
i SA TO ya 
Bet CALL MIYTAP ¢x,Y) 
2 BR CANTINVE 
Ba t CA TINUE 
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a a TES eS ee ee ‘oy SETS ARE 


34 


364 


369 
370 


368 
365 


375 


376 
377 


oF 


465 


PT us 


a eae 


STATE 


Cou 
C4utr 
Cau 
Cate 
Cau 
CALL 
C&L 
DO 34 
i ae 
x 3s ¢ 
yso0 
CALL 
x = YX 
Yay 
CALL 
CANT] 
KT = 
Ge 36 
KT = 
IF (« 
Fac s 
KT s 
69 19 
FAC 
R23 
x sR 
IF (xX 
Ys 0 
CALL 
Yee 
CALL 
CSNTl 
KS s 
KJ!” = 
OS 37 
«§5 = 
EF CK 
Fac = 
«S s 
G3 TS 
FAC s 
Hos & 
Yosr 
x £2 0 
CALL 
a a) 
CALL 
ceuty 
C4uL 
IF (% 
SFACS 


MENT 


STICK (1,26,5,5,5) 
ATIC COT FloteSetO?) 
ATIC«K(73,141,2,%9) 
ATICK (1 S2.151 et %e5) 
OTICK(1S6,146,5,",2) 
ATICK(168,178,251,2) 
ATIC# (179,181 ,1e1e2) 
kKF Ss ly 6 
JACKE) 
oD 
et 
PLEeT(X.Y,3) 
X ONE) 
Y (NF) 
PLOT(CX,Y22) 
NUE 
9 
uo«Y 3 {, 400 
Kr + { 
Tome. to) Go Te sas 
209 
0 
379 
1,9 
KY © | 
SexMAKSEMAX 
oT e CXMAX @ 69991)) GA TS $65 
Cau 
PLOT(X,Ys3) 
Fac # DEL 
PLST(xk,%,2) 
NUE 
9 
HMAX / SINT © 1,004 
KJ sl, kM 
4S + | 
S.%E,10) GP TS 374 
2.0 
a 
377 
1.9 
INT « (KJ 01) 
RYAAXSHMAY 
ec 
PLOT (X+Y,3) 
FAC * DEL 
PLOT(CK,Y,2) 
SUE 
PENCRKG(11) 
MAX@eY”AX) 460,460,461 
xMAX 4125 
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a a STS ee Smee eee. er nn Make 


“ee 


sail 


= eS ne a ee ee 
7 


TOSS 


Raina 


s 
q 


So! 


402 


Bul 


a3 


STATEMENT 


G* T9 &Bae 

SFacsy“ske 125 

re,i75eSF 4c 

o9 100 [* 3 1, 5 

Nos Ri Ac ¢ TREAC (IS) 

IF (NRA GT.I0) GS TS 104 

TRUNTT s&s TEFaAC (TR) 

69 Te 102 

CAXTIRVE 

TRENT TSTIFPESC(S) 

RG to £4 2 fe & 

‘im So mMAN / TREAC (Ih) 

IF (NH,GT.10) GO TS 110 

TeuSTT = yer ac (Ih) 

a3 Te 128 

CArrInve 

THUMITSImF Oc (8) 

ase, 05eSF aC 

You ,SeSFac 

CALL TNUMPM (X,Y,4,9,0,0) 

" = TRENIT 
IDGITSS4L OG OCFLOAT(N) )41,00000! 
C4LL SENTER (m,ICGITS,IOGITS,874S) 
X & (N/R“AX) w@ XMAX @ KITAS 

JF €X% + SIAS ,GT,. X%AX) GO Te EOD! 
CALL TNUMBE (X,Y —h pn, 0,0) 

“ SN ¢ [TRUNIT 

68 T9 129 

Yset,o#SF ac 

CALL SENTES CH, 21,21,8148) 

X¥ 5s U,8 # XMAX @ STAS 

CALL LETTER (X,Y,m,2thRANGE, “NAUTICAL MILES,9,°,21) 
Xse,SeSFAC 

Ysa, 9875 *SF AC 

CaLh INUMRE (X,¥,h,0,0,)) 

» S Teun] T 
THGITSEalL Shi c CFL SAT (ND) 43, 000001 
KeeC IS#e#ITHGITSa#SFAC 

YEN /SHMAX)HYMAKO,UATSeSFAC 

Tr ¢Y¥ ‘ te VRE} G* TS 803 

C4ubL INUMRE (X,Y,H,N- 0.0) 

* 3s N # Jmun{[t 

Ge 79 (21 

RS ef ,40eSFAC 

CALL SENTER(*,33,33,F14S) 

Ys 0,5 * Y*&X © 94s 

CALL LETTES (X,Y ph, 33R TARGET HEIGHT ABSVE ANTENNA, FEET,90,,33) 
x{Sxe 3, er 

CALL SENTER (HK, 26,26,R14S) 

CALL LETTEE (X1,% 4, 2en(4NTENNS HEIGHTS FEET) ,90,,26) 
Y2SYe1U,S7 {en 
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em wt 
Taki ineitt Me 


aay 


aa7 


STATEMENT 


CALL NIMSER (K1,Y2,4,A4NTHGT,90,,01) 
XPQS 2 ,exdn 

YP sed yen 

ce 6ae IL & t, 2 

INO & TANGCIL) 

Aass4aa.GCIL) 

Cxs,1 

TF CIL.GE,6)CXS,125 

XEXKCISO) © .GeCTICINO) eSFAC @ CxwSFAC 
YEYYCIND) & g4eSNLCIND #SF4C © ,O87SeSFAC 
IF (IL ,€G, 9) X 2x eb 


IFCIL,ED,% o9R. (YMYER) GTQ(1,5e4)) GE TH ERG 
IF (CKPReX) LT, (CU, 0eF) 8h, KL T,2,e%) GM TS BNL 


xPUsx 

ypRsy 

IF CYL.GE,e) GS TO 888 

CALL NUMFES (X,Y¥,m,44%,0,0e1) 
xs X © ,15 # SFAC 

6° TS asg 

CALL NUMBES (%,Y,%,4A%.0,9-01) 
CeMTINUe 


CALL DFGSER (KX @ ,35eSFAC, Y + 4175 # SFAC,,OA#SFAC) 
3 CONTINUE 


END 
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oy ae perros = <= Sib a aT idk ee TR 


om i ; 


piers 


Se eileen, en he coi 


STATEMENT 


FusOTyen Fic) 

CIM AR sdf HEF, GRAD, FAD, CANST, 10182), 
yA SREP RE KP (GRAN eK) 

COSXs24n 

VS2,0 RR eERRSR 

as2,uNlCeoCl eC 

Fx 3 SORT (U(N) © V @ * @ Ven) 

et. 2 CRNST &@ C1,0 eVIRCL.0 + CO)/EX 

ENC 


39 


N 


STATEMENT 


SUBRSUTINE SIMCA (KL, XEND, TEST, LIM, SRE4,NS1,%,F) 
apace 
foNso,4 
INTS1 
; Vst,.0 
F EVENSU,9 
; AREAL30,0 
19 ENDSSF (KL) OF CKEND) 
2 MS(KXEND eX1L)/V 
S0DSEVE'+8CD 


XBX1+H/2, 
i EVENS0,9 
, OO 3S Tst,1NT 
a 2} FVEN SEVENSF (Xx) 
x ABXor 


3 CONTINUE 
* 31 AREAS(ENDSO4 OWE VENS2 04200) ahsO,f 


NOTSNOT OL 

: 34a RBABS (CARE ALCAREA)/4REA) 

4 TF (NS ToL IM) 341,35,35 

r. 341 IF CReTEST)35,35,4 

ie 3s RETUPK 

3 4 AREALZAREA 

bn 46 INTS2ea INT 

‘4 Va2,0"V 

th Ge Te 2 

3 END 
* 

2 

ii 

3 

¥ 

f 

t 

5 

¥. 

# 

a 

40 

ech | 

«: 

— ee ea ee a os ‘oy, GEO OL SSS 
mien tae . oe ont — ref 2 gE ie ie Pa a i st aii ie 4 + her . . : 


: ——— samc 


STATEMENT 


SrRROUTTE ATIC (Ta, JS, Ka,MF,MC) 
Cees ses eeltAry, YYCLAL), CTLCLEN), ShrciAr), CEL 
“a s *E 
Fe jm 3 Th, Ja,Ks 
MA sz “A @ fj 
TP C48 NEMEC) GA Te S 
2 Fac = 2,0 
MA zs” 
Ge Tf 
FAC 3 3,0 
ao XK S Rx) 
Ye V¥C(K) 
Cane PL ET (X,Y 35) 
x eX @ CEL eFACRCTICK) 
Ys Y » CEL aFAC#SNI (kK) 
Cath PLAT(X,Y22) 
1 CentTInve 
ENC 


et 


ore 


STATEMENT 


SUBRPUTINE SENTER (h,4,2,824S) 
CwenST s ,28S7143 2 
ato Tae (3ele}) « CONST 


aIDTeA S (Berel) # C8NST 
ATOTHE © 0,5 # wICTHN 


Fr4s s 
RE Tm 
BRD 


STATEMENT 


Su@@9uTI*e wINTAe (X,¥) 

CAMMAN SMTFES X25 V2, Ke Vip KA, YO, ERA 
COMMON SEAS MQ 

Cara sos 

Oata usos 

OXCUL,Vi,Le,V2) 2 SGRTC('aeuUl)ew2 + (V2—eV1) eee) 
O1 3 Ov€xt.¥1,x%2,¥2) 

TF (09 ,E9, 9.) GY TA II 

C2 8 Cx(x,¥,X1,¥1) 

TF (02 EG, 0.) GY TC 2 

C3 = Ox(x,¥,x2,¥2) 

Pty S&S CX(xb,VA, 2,72) 

IF CO} obF. CLV) GO TO 4 

CISN SF COSCO V*eHSel 1 eO1 js (2, ef 1902) 

IF (CeS% ,GT, 1, 3%, COSN LT, of.) CASA = 1, 
STINE &S SGRT(1, © CI§N#C%SN) 

CEVN & D2eSIté 

TF (COEVN LE. ERROR) G* Te 2 

COLL PLST (X4,%a,2) 

es Mel 

x2 mA 

v2 Ya 

x} 

vi 

xa 

ya 


ne Seager Raa eae, of pa a ee 
ORE Cre tel ree 2 My ARMM RRO i ti ee Atta te Gin 4 . 


STATEMENT 


SLESELTINE SRAREF (6%HZ, PST, YPEL, Sho, PHI) 
CeMPLEn EPSC, GA%, SGTERM, TER» 
Darts PLAST 4 O,/ 
SINPSI = SIN (PST) 
CSPSI = COS(PSI) ee? 
IF (FMHWZ Eu, FLAST) G® TR 200 
FLaST s FMnHZ 
» & 299,793 ss FMR] 
IF C(FMNZ GT, 1500,) GS Te 151 
\ 1$0 SIG 2 4,3 
| E°Ssi = 3G, 
ce Te 455 
151 SIG © 46.3 ¢ (FMHZ © 1500,) * ,00148 
IF CFMRZ .GT, 3090.) GA TS 154 
153 EFS! s 80, © (FMH2 © 1500,) * ,00733 
ce Te 1558 
1546 EPS! 3 69, e(FMHZ © 3040,) * 6002429 
SIG 2 6.52 * (FMHZ @© 3000,) * ,001314 
E - {55 EPSC 3 C¥PLX (EFPS1,"60,"4wSIG) 
F 200 SGTESM= CSRRTCEPSCeCSPSI) 
IF CIPRAL oNE, 1 ) Go Tf lol 
169 TERM s EPSC # SINPST 
f GAM & (TERM@SGTERM)/(TERMSSQOTERM) 
¢ Ge TS 180 
161 GAN & (STNEST @ SUTERM) / (CSINOST + SGTERM) 
180 Gre =s CaRS (Gam) 
Pet SedTAn? (ATMAG(G4M), REAL (GAM)) 
PETURA 
K ENG 
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. a — 
aeeme- MEET - ; 


Al Pn. he Pee Le eo eS 


ioe ree es a en ‘ ‘yg, RNS i eS 


STATEMENT 


SUBROUTINE OFGREE(X,Y,e) 
MDIMENSTON CELK(6), OELY(3) 
DATS DELO TOT pele Ope TOTL Ve O ee TOMI ol oO He Me 707I 49,07 
DATA DELYSO TOT Ce Dp eO TOTS pel al eG gM ep Delsey 1y,0/ 
hos Re 4ju2is 
ks Ko Sen 
Ys Ve ,5#0 
CALL PLOT (X,Y,3) 
Of 1 Ir si, 8 
KX =X # OELXCII) eb 
Y 2 Y + DELY(II)#D 
{ CALL PLAT(X,Y,2) 
END 


Bake TOR 


Bs 


aie 


! 


45 


re ee eee anienaka . - 1 ge Os 2 ARETE 
em ew eS a Ee ro ‘ woe otealemnementl ae “te? a , 


rae RABID NR ERT TE I Ie ag MES Ee ARNE PRR ne CME I I Bree 


STATEMENT 


SUBRECTINE ISTRST CxL,¥1,X2,V2,Xh,¥b,XB,YB,x0,0) 
TF (ABS(xZext), GE,10,be76) GY TO 2 
{ IF (ARSC keeKa) LT, Lu,E 76) GS TH 99 
sO 3 KY 
S$ s (YhevA)/(XRexa) 
YO 3 Sw(¥Qex4) + YA 
BE TIRN 
TF (45S (xGexk) GE,10,EF"76) GS TH 4 
IF CARS (XK PexKt) LT,10,E976) G9 719 99 
»O S&S x4 
$ = (v¥20ev{)/(x2e%}) 
YQ & Sef ¥Cex{) + YI 
rE TUEN 
a St Ss (Y¥PeVA)/S(XBexKA) 
$2 = (v2eVt)/(K2ex1) 
IF ¢St «EG. Og) Gh Te § 
IF CABSCS2@S1) LT.10,E976) GF 18 99 
eT® = 52/81 
YO © (v2eSew(xbeox2)evVsakTG)/(},e8T8) 
xO 3 (Y¥IeV4)791 6 XA 
QE TURN 
§ IF (S82 «EQ, Oey GO TS 99 
¥Os (YAeY¥2)/S2 + x2 
YOSYA 
RETURN 
Q9 FRINT 19°C, X1y YI 
PRINT Jul, ¥2,7%2,%4,Y4,X8,Y8 
x0 2 0, 
vo 3 0, 
100 FUBMAT (3¥,O4RKCALL TE ITNTESECT ARGFTED, LINES PAPALLEL, NO INTERSE 
wOTION, XE FS GED eeyerm YE S GELS 42) 


nN fh 


4 101 FARMAT( SX, 6H X2 SF pEL0,2,6H Y2 = E1992, 
# OM KA FS -E10,2,6h YA = 4£10,2,6% XB F 4E10,2,6" YR = 4E10,2 //) 
END 


i 
: : 
: 


STATEMENT 


SURREUTINE APRS (X1,%1,x2,%2,TIPL) 
CAYMAN SAG/SANG 

OaTa ANG / ,35/ 

YA z= ye @ Yj 

x4 3 x2 © KL 

A s ATAN2(YA,X4) 


43 2 A * ANG 
44 3 A @ ANG 
X33 x2 © TIPL * C9S(43) 
¥3 = ¥2 © TIPL ® SIN(43) 
XG 3 X2 © TIPL #* C9S( AG) 
Yu s y2 © TIPL * SIN(44) 


CALL PLOT (Xi, Yh, 3) 
CALL PLAT (xe, Yeo 2) 
CALL PLOT (x3. Y3, 2) 
CALL PLAT (x4, Yu, 3) 
CALL PLOT (x2, Y2, 2) 
END 


STATEME'T 


SULERSUTINE REZERS 
ENT@Y PEN. CHG 
i ENTRY DASHAN 
ENTRY DAaSn4E 
END 


; 
: 
© 


= 


APPENDIX B: Input Data Card Formats for LOBMUF (For ASC - Jan 1977) 


Program LOBMUF will plot the envelope of N radars using the pattern 
calculations contained in the original program described in NRL Report 
7098 by L. Blake. Incoherent integration is assumed. An antenna height 
must be specified in subroutine RHACHT. The parameter name for this 
antenna height is ANTHITE and is used to draw the range-height grid. The 
actual antenna heights for the individual radars are input separately and 
used for the multipath calculations. ANTHITE isnot very important except 
for the labeling of the height axis of the grid. 


Input Card Sequence 


Col. 

1. ANTHITE in F10.0 field 

2. Label Card 

3. Grid Parameters in F10.0 fields: 
a. XMAX 1-10 
b. YMAX 11-20 
c. RMAX 21-30 
d. HMAX 31-40 
e. MIN 41-50 
f. THMAX 51-60 
g. WHFT 61-70 
h. RDR 71-80 


= maximum X-dimension of chart in inches 
YMAX = maximum Y-dimension of chart in inches 
RMAX = maximum range of chart in nautical miles 

= maximum height on chart in feet 

THMIN = minimum elevation angle for the plot 

THMAX = maximum elevation angle for the plot 

WHFT = sea wave height in feet 

RDR = number of radars for which an envelope is to 


be plotted. 
4. Common Parameters in F10.0 fields: 

a. PDT 1-10 

b. PFA 11-20 
c. PULS (total number of pulses) 21-30 
d. CASE 31-40 
e. AHFT 41-50 
£. TILT 51-60 
g. POL 61-70 
h. CSC 71-80 


PDT = probability of detection 

PFA = false alarm exponent, i.e., the positive value of 
the exponent (power of ten); for 107°, enter the 
number 6.0, etc. 

PULS = total number of pulses 


48 


Se RS ae 


eo 


CASE 
AHFT 
TILT 


Swerling case number 

antenna height in feet 

tilt angle of the antenna beam maximum with 
respect to the horizon in degrees. 


iow 


For reference, POL is polarization as follows: 


POL ay, vertical 
POL = 2, horizontal 


CSC is for pencil or cosecant squared antenna pattern: 


csc = 0. pencil beam 
csc = 1. cosecant squared beam 


5. Radar Parameters in F10.0 fields: 


a. RFS 1-10 
b. FMHZ 11-20 
c. BWD 21-30 
d. SLDB 31-40 
e. PULNUM (number pulses at FMHZ) 41-50 
£. FREF 51-60 


Repeat card 5 for N repetitions. 


RFS = calculated or assumed free-space range of the 
radar on the specified target 

FMHZ = radar frequency in megahertz 

BWD = antenna half=-power beamwidth in degrees 

SLDB = first elevation sidelobe level relative to the 
main lobe 

PULNUM = number of pulses at each specified frequency 

FREF = frequency used to calculate RFS if RFS was 
calculated for the total number of pulses (PULS) 

FREF = 0, if RFS was calculated at each frequency using 
PULNUM. 


After cards 1-5, the sequence can be repeated as many times as desired, 
starting with card 2. 


A note of caution is that the elevation plot angle increment is 
determined by 


THMAX ~ THMIN 
2000 


The value has proven adequate for S-band radars where 10° and 0° were the 


values for THMAX and THMIN. It is recommended that an envelope for a 


Single radar (RDR = 1) be run when there is doubt that this angle increment 
is small enough. 
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